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Abstract 



We investigate the dynamical stability of a fully-coupled system of N inertial rota- 
tors, the so-called Hamiltonian Mean Field model. In the limit iV — > oo, and after 
proper scaling of the interactions, the /i-space dynamics is governed by a Vlasov 
equation. We apply a nonlinear stability test to (i) a selected set of spatially ho- 
mogeneous solutions of Vlasov equation, qualitatively similar to those observed in 
the quasi-stationary states arising from fully magnetized initial conditions, and (ii) 
numerical coarse-grained distributions of the finite- iV dynamics. Our results are con- 
sistent with previous numerical evidence of the disappearance of the homogenous 
quasi-stationary family below a certain energy. 
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1 Introduction 

Nonextensive Statistical Mechanics (NSM) was proposed in 1988 [1] as an 
attempt to account for many phenomena that can not be described within the 
standard Boltzmann-Gibbs scenario (see [2] for reviews). NSM starts from a 
generalized entropy functional, and using suitable constraints, a generalization 
of the usual canonical ensemble thermostatistics is obtained. Instead of the 
exponential weight, in NSM one has a g-exponential probability distribution 
in T-space: 



(if the quantity between curly brackets is negative then P = 0). Here, x is a 
point in the iV-body phase space of a system described by the Hamiltonian 
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P(x) oc expj-/3#(x)] = {1 - /3(1 - q)H(x)} 



1/(1-9) 
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H(x), (3 a generalized inverse temperature and q the entropic index [2]. Unlike 
the canonical distribution function [Eq. (1) with q = 1], which applies to 
equilibrium situations, the g-canonical distribution (with q ^ 1) is believed to 
be associated to meta-equilibrium or even out-of-equilibrium regimes [2]. 

The predictions of the theory can in principle be derived from Eq. (1). How- 
ever, explicit calculations for the anomalous high-dimensional systems of in- 
terest for NSM have not been performed yet due to large technical difficulties. 
For instance, single particle distribution functions require the integration of 
P(x) over N — l particles. 

On the other side, there are many experimental and numerical data well de- 
scribed by g-exponential functions [2]. This agreement is considered as an 
indirect evidence of the applicability of the NSM to those systems. More- 
over, g-exponentials have also been observed in non-Hamiltonian and even 
non-physical (biological, economical, etc.) systems [3], but in these cases the 
connection with Statistical Mechanics (from first principles) is less clear. 

A model system potentially important for NSM is the so-called Hamiltonian 
Mean Field model (HMF), an inertial XY ferromagnet with infinite-range 
interactions [4,5]: 

i N I N 

where the coordinates of each particle are the angle 0i and its angular mo- 
mentum pi. This model can be considered a simple prototype for complex, 
long-range systems like galaxies and plasmas. In a certain sense, the HMF is 
a descendant of the mass-sheet gravitational model [5]. 

The HMF has been extensively studied in the last few years (see [6] for a re- 
view), especially due to its long-lived quasi-equilibrium states, which exhibit 
breakdown of ergodicity, anomalous diffusion, aging, and non-Maxwell velocity 
distributions. In some cases g-exponential behaviors were reported: two-time 
correlation functions [7,8], single particle velocity distributions (truncated q- 
Gaussians) [9] and temperature relaxation [10]. The simplicity of the HMF, 
that made possible a full analysis of its statistical properties both in the stan- 
dard canonical [5] and microcanonical ensembles [11], gives support to the 
hope that it may be possible to derive explicit first principle predictions, e.g., 
for single particle velocity distributions. 

The most widely investigated quasi- stationary states (QSS) of the HMF arise 
from special initial conditions (water-bag) that evolve to states with well de- 
fined macroscopic characteristics and whose lifetimes increase with the system 
size [9]. For a given energy, the QSS is different from the corresponding equi- 
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librium state to which the system eventually evolves. These states have been 
observed for specific energies e = H/N below the critical value e CT = 0.75 cor- 
responding to a ferromagnetic phase transition. In a magnetization diagram 
(see [9] and Fig. 1), a first part of the QSS family lies on the continuation 
of the high-energy equilibrium line (having zero magnetization) down to the 
energy e min ~ 0.68. At e min the m = line ceases to exist, being substituted 
by a curve of magnetized states. 




Fig. 1. Squared magnetization versus energy for quasi-stationary states (circles). 
System size is N = 10 4 . We averaged over 10 initial conditions and in the time 
window 500 < t < 2000 (see Sec. 4 for numerical details). The full line is the 
analytical result for thermal equilibrium [5]. 

This paper is an attempt to understand the origin of such discontinuity. We 
base our analysis on the associated Vlasov equation which describes the contin- 
uum limit of the Hamiltonian particle dynamics. Using the stability criterion 
recently developed by Yamaguchi et al. [14], we show that the coarse-grained 
homogeneous (zero-magnetization) QSS family becomes unstable at e m i n . 



2 Vlasov approach 

Introducing radial and tangent unitary vectors = (sin 6>j, cos 6i) and 0* 
(cos6>j, — sin#j), Hamiltonian (2) leads to the equations of motion 



0i=P» 

Pi = m-6i, (3) 
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for i = 1, . . . , N, with the magnetization per particle 



m = — 




(4) 



Eqs. (3) describe the two-dimensional motion of a cloud of N points in the fluc- 
tuating mean-field m(t). In the continuum approximation, the cloud becomes 
a fluid governed by the Liouville equation for the single-particle probability 
distribution f(6, p, t) 



Liouville equation (5) has to be solved self-consistently with the calculation of 
m(£), and is thus formally equivalent to the Vlasov-Poisson system found in 
plasma and gravitational dynamics [15]. The continuum approximation will be 
progressively better as N grows, especially when the cloud of particles is spread 
over a finite region in /x-space. This is the case of the distributions we will 
analyze in this paper (see Fig. 2). Detailed discussions of Vlasov approximation 
can be found in [12,13,14]. 



3 Dynamical stability 

Every spatially homogeneous distribution f{p) is a stationary solution of 
Vlasov equation (5). The reason is that homogeneity leads to m = and thus 
momenta do not change. The interesting point is whether these homogeneous 
solutions are stable or not. There is the well known Landau analysis [15] which 
concerns linear stability. This linear criterion was recently applied to the HMF 
by Choi and Choi to discuss ensemble equivalence and quasi-stationarity [16]. 
Besides the linearity restriction, Landau analysis has the disadvantage of re- 
quiring the knowledge of detailed analytical properties of f(p). This excludes, 
in principle, numerically obtained momentum distributions. 

A more powerful stability criterion for homogeneous equilibria has been pro- 
posed by Yamaguchi et al. [14]. This is a nonlinear criterion specific to the 




(5) 



with 




(6) 
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HMF. It states that f(p) is stable if and only if the quantity 



— oo 



is positive (it is assumed that / is an even function of p). Note that this 
condition is equivalent to the zero frequency case of Landau's recipe [5,16]. 

The authors of Ref. [14] showed that a distribution which is spatially homo- 
geneous and Gaussian in momentum becomes unstable below the transition 
energy e cr = 3/4 (see also [4,16]), in agreement with analytical and numeri- 
cal results for finite N systems. They also showed that homogeneous states 
with zero-mean uniform f(p) (water-bag) are stable above e = 7/12 = 0.58... 
(see also [5,16]), and verified that the stability of these Vlasov solutions has a 
counterpart in the associated particle dynamics. 

In the same spirit, it is also worth analyzing the stability of the family of 
g-exponentials 

f(p) oc exp g (-ap 2 ) = 



1 - a(l - q)p> 



(8) 



which allow to scan a wide spectrum of probability distributions, from finite- 
support to power-law tailed ones, containing as particular cases the Gaussian 
(q = 1) and the water-bag (q = — oo). In Eq. (8), the normalization constant 
has been omitted and the parameter a > is related to the second moment 
(p 2 ), which only exists for q < 5/3. In the homogeneous states of the HMF one 
has (p 2 ) — 2e — 1, as can be easily derived from Eq. (2). Then, the stability 
indicator / as a function of the energy for the g-exponential family reads 

* = 1 - "7 ^TT V • (9) 

2(5-3g)(2e- 1) V ; 



Stability occurs for energies above e 



q 



e q (q) = e CT + ^f-^ . (10) 



It is easy to verify that one recovers the known stability thresholds for the 
uniform and Gaussian distributions, i.e., e q (— oo) = 7/12 and £ q (l) = 3/4. 
We remark that Eq. (10) states that only finite-support distributions, corre- 
sponding to q < 1, are stable below e CT . This agrees with numerical studies in 
the QSS regimes of the HMF (see Fig. 2). 
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Fig. 2. Momentum distributions for energies indicated on the figure. Histograms 
(bin size ~ 0.07) were generated averaging on the time interval (500,2000) and 10 
different initial conditions (very small random perturbations of the fully magnetized 
regular water-bag). In all cases N = 10 4 . See Sec. 4 for further details. Full lines are 
water-bag-plus-cosine fittings to the data. 

The QSS we discuss in this paper are obtained from particular initial condi- 
tions: all angles are equal and momenta distributed uniformly with zero mean; 
energy is chosen slightly below e CT . After an initial "violent relaxation" (the 
initial state is highly inhomogeneous) the system evolves to a QSS charac- 
terized by spatial homogeneity and momentum distributions with cutoffs [9]. 
However, the resulting distributions of momenta are not uniform any more. 
Some structure appears on top of a slightly smoothed water-bag (see Fig. 2). 
A very simple family of functions exhibiting the basic features of the observed 
f(p) is the water-bag plus cosine: 

ft \ _ / V( 2 Pmax) + -Bcos(7rp/p max ), for 

0, otherwise. 1 j 



The stable functions f(p) occupy a region in parameter space (p max , B) de- 
fined by three boundary curves (see Fig. 3). The left-top boundary is given 
by I = and the other two are associated to positivity, i.e., f(p) = 0. The 
metastable states obtained from numerical simulations of finite systems (to be 
described below) can be approximately fitted by water-bag-plus-cosine func- 
tions. Remarkably, the resulting pairs (p max , B) fall close to the boundary of 
Vlasov stability, and exit the stability region for energies below the limiting 
value e ~ 0.68. 

Fig. 3 also informs that for low energies, e < 7/12, there is a qualitative change 
in the shape of stable homogeneous distributions: they change the maximum 
at p = by a minimum since B < 0. As e — > 1/2, the stable distribution 
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Fig. 3. Stability diagram for water-bag-plus-cosine momentum distributions. 
The shadowed area corresponds to the stability domain. Also shown are three 
iso-energetic subfamilies with e = 0.50, 0.60, 0.69 (dotted lines). Symbols correspond 
to quasi-stationary states obtained numerically, for energies e = 0.67, 0.69, 0.72 and 
N = 10 4 (see Sec. 4 for numerical details). 

resembles a couple of delta functions approaching the origin. 



4 Numerical results 



The application of the stability criterion to a discrete distribution requires 
some coarse-graining. We will consider smoothed distributions 



1 



A? 



iv i=l 



(12) 



with Lorentzian delta functions 



5 1 (x) 



7 1 



7T X 2 + 7' 



2 ' 



(13) 



The coarse-graining parameter 7 corresponds to an intermediate scale, be- 
tween microscopic and macroscopic (to be specified). We will also impose 
parity, f(p) = f(—p), so that our smoothed distribution will be calculated as 



A(p) = ^ E 



+ 



(P-Pi) 2 + 1 2 {.P + Pif + I 



(14) 
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Accordingly, the stability indicator becomes simply 

^ 1 + fe(frw' (15) 



All simulations have been done with N = 10 4 particles. We have considered 
water-bag initial conditions. All particles are set at 6 = 0, corresponding to 
maximum magnetization modulus m = 1. Momenta were distributed on a reg- 
ular lattice, symmetrically around p = 0. Our experience indicates that this 
"crystalline" distribution behaves like a larger system, say N = 10 5 , where 
the p's are extracted randomly from a uniform distribution. Initial states were 
propagated according to Hamilton equations, approximated by a symplectic 
fourth-order scheme [17]. We registered the quantities of interest m(t) (homo- 
geneity indicator) andpj(t) (input to the stability indicator). Strictly speaking, 
m = does not imply that the states are homogeneous. However, the sudden 
relaxation that leads to the present QSS mixes particles [8] in such a way that 
m = and spatial homogeneity are expected to be synonymous. 

Concerning the calculation of smoothed distributions, we chose a coarse-graining 
parameter 7 = 0.05, after verifying that such a choice erases microscopic 
fluctuations, but preserves macroscopic features of the QSS momentum dis- 
tributions (like those shown in Fig. 2). Anyway, some tests were run with 
7 = 0.025 and 7 = 0.1 to check that our results are not substantially changed 
by the precise choice of the smoothing (see inset in Fig. 5). We also eliminated 
rapid time fluctuations by doing an additional average in the time window 
500 <t< 2000, a typical domain of quasi-stationarity after initial relaxation 
of the water-bag (see Fig. 4). 

The sensitivity to initial conditions (of the microscopic state of the system) 
after times like those considered here (t « 1000) is enormous. To be sure 
that the effect we are observing is statistically robust, we also did an average 
on initial conditions. These were obtained by slightly perturbing the regular 
water-bag with a small random component: Pi = Pi(l + 10~ 4 £), where £ is 
random uniform in [—1,1]. Our stability analysis of the coarse-grained quasi- 
stationary states of the HMF is summarized in Fig. 5. There one sees that 
stability is positive for energies above e ~ 0.68, except for the small region 
0.73 < e < 0.75. We believe that this is probably a spurious effect due to the 
proximity of the ferromagnetic phase transition; further analyses, involving 
longer times, or larger system sizes, are necessary for elucidating this ques- 
tion. The fact that the stability indicator becomes negative below e ~ 0.68 
signals the disappearance of the homogeneous metastable phase at that en- 
ergy. Simultaneously with / crossing zero at that point, the magnetization 
rises from negligible values to finite ones, of the order of those corresponding 
to equilibrium (Fig. 1). This stability picture is consistent with the analysis 
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Fig. 4. Example of stability as a function of time, for the fully magnetized regular 
water-bag initial condition, with e = 0.69, N = 10 4 , and smoothing 7 = 0.05. The 
full line corresponds to a running average over a time window of width At = 100, 
indicating quasi-stationarity in the interval 500 < t < 2000. 

of the water-bag-plus-cosine family (Sec. 3), showing that distributions like 
those observed in QSS live close to a stability border, and that this border 
is traversed when energy goes down a limiting value. The symbols in Fig. 3 
were obtained by fitting a water-bag-plus-cosine to the numerical histograms 
of Fig. 2. 

A comment about negative stability is in order. The present stability test 
only applies to homogeneous states. However, below e = 0.68 the measured 
distributions are evidently inhomogeneous (m 7^ 0). In these cases, negative 
stability refers to hypothetical homogeneous states having the measured f(p). 



5 Conclusions 

Yamaguchi et al. conjectured that stable homogeneous solutions of Vlasov 
equation lead, via discretization, to quasi-stationary states of Hamiltonian 
dynamics, and showed that it is indeed the case for some families of states of 
the HMF [14]. We have shown that the reciprocal may be also true: metastable 
states of the discrete Hamiltonian dynamics lead, after smoothing, to stable 
solutions of the Vlasov equation (this was verified for the homogeneous QSS). 

Moreover, our analysis indicates that the disappearance of the homogenous 
family (resulting from the violent relaxation of special water-bag initial con- 
ditions) at e ~ 0.68 can be interpreted as a reflection of the lose of Vlasov 
stability at that critical energy. This interpretation is supported by similar ev- 
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Fig. 5. Stability versus energy for quasi-stationary states. Note that I goes through 
zero around e ~ 0.68. The smoothing parameter is 7 = 0.05 (N = 10 4 ). We averaged 
over 10 initial conditions and in the time interval 500 < t < 2000. Inset: I vs e for 
7 = 0.025 (V), 0.05 (•) and 0.1 (A). 



idence resulting from the analytical study of the water-bag-plus-cosine family. 

Despite providing interesting information, the results we presented are not 
definitive. Some points remain unclear and deserve further investigation, e.g., 
the stability properties close to the ferromagnetic transition. Another point, 
perhaps more important, concerns the quantification of spatial homogeneity. 
A null squared magnetization indicates the absence of spatial structures on the 
largest scale. A more stringent test would require the verification that higher- 
order Fourier coefficients of the particle density tend to zero in the large N 
limit. 
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